por Francisko de Moraes Rezende, , 2020-03-26

# carregando os pacotes ---------------------------------------------------

lista_de_pacotes  <- c("tidyverse",
                       "readxl",
                       "here",
                       "lme4",
                       "janitor",
                       "ggbeeswarm",
                       "broom",
                       "tufte",
                       "ggthemes",
                       "ggpol")


pacotes_novos <-
  lista_de_pacotes[!(lista_de_pacotes %in% installed.packages()[, "Package"])]
if (length(pacotes_novos))
  install.packages(pacotes_novos)


library(tidyverse)
library(readxl)
library(here)
library(lme4)
library(janitor)
library(ggbeeswarm)
library(broom)
library(ggthemes)
library(ggpol)

Carregando os dados

ANOVA ui

Código

ui_m1 <- lm(ui ~ sp * face, data = ui)

anova(ui_m1)

ui_m2 <- update(ui_m1, . ~ . -sp:face)

anova(ui_m2) %>% 
  rename(`P value` = `Pr(>F)`)

Tabela ANOVA ui

anova(ui_m2) %>% 
  rename(`P value` = `Pr(>F)`) %>% 
  knitr::kable(digits = 3)
Df Sum Sq Mean Sq F value P value
sp 1 1.082 1.082 238.52 0.000
face 1 0.000 0.000 0.07 0.798
Residuals 9 0.041 0.005 NA NA

Gráfico

Os dois gráficos abaixo tem as mesmas informções, escolha o que mais lhe agradar.

ui %>% 
  ggplot(aes(x = sp, y = ui)) +
  geom_boxplot() +
  geom_quasirandom(groupOnX = T,
                   alpha = .5) +
  theme_classic() +
  theme(axis.line.x = element_blank(),
        text = element_text(size = 12)) 

ui %>% 
  ggplot(aes(x = sp, y = ui)) +
  geom_boxjitter(jitter.height = 0) +
  theme_classic() +
  theme(axis.line.x = element_blank(),
        text = element_text(size = 12)) 

ANOVA angle

Código

angle_m1 <- lm(angle ~ sp*face, data = angle)

anova(angle_m1)

angle_m2 <- update(angle_m1, . ~ . - sp:face)

anova(angle_m1, angle_m2)

anova(angle_m2)

Tabela

Df Sum Sq Mean Sq F value P value
sp 1 11803.277 11803.28 586.646 0.000
face 1 17.280 17.28 0.859 0.378
Residuals 9 181.079 20.12 NA NA

Gráfico

Os dois gráficos abaixo tem as mesmas informções, escolha o que mais lhe agradar.

angle %>% 
  ggplot(aes(x = sp, y = angle)) +
  geom_boxplot() +
  geom_quasirandom(groupOnX = T,
                   alpha = .5) +
  theme_classic() +
  theme(axis.line.x = element_blank(),
        text = element_text(size = 12)) 

angle %>% 
  ggplot(aes(x = sp, y = angle)) +
  geom_boxjitter(jitter.height = 0) +
  theme_classic() +
  theme(axis.line.x = element_blank(),
        text = element_text(size = 12))